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Abstract 



We show that quantum computation can be performed in a system at 
thermal equihbrium if a spontaneous symmetry breaking occurs. The com- 
puting process is associated to the time evohition of the statistical average of 
the qubit coherence operator. This average defines logical states which evolve 
under the action of impulsive gate perturbations. Non trivial quantum coher- 
ence properties appear in the ordered phase characterized by a spontaneous 
symmetry breaking. Qubits are associated to spins and interaction with the 
lattice vibrations is approximated by an effective spin-spin interaction. A 
finite decoherence time is obtained only if spin correlations are taken into 
account in the framework of the Ising model as corrections to the Mean Field 
approximation. Decoherence is slowed down as the order parameter saturates. 
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In the quantum computing process information is stored in qubits i.e. in a quantum 
state of a two levels system. Coupling with the environment produces an "entanglment" 
of the qubits with the degrees of freedom of the environment. This "entanglement" is the 
origin of a progressive loss of information, a process known as decoherence. Hence, systems 
at thermal equilibrium have not been considered suitable for a computing process. We show 
that quantum computation is however possible even in thermal equilibrium in presence of a 
spontaneous symmetry breaking. Logical states of computer will be associated to suitable 
thermal averages. The 0-logical state is associated to thermal equilibrium and 1-logical state 
to local excitation determined by an impulsive perturbation. Decoherence is then associated 
to the local excitation relaxation toward equilibrium. We shall see on a specific model 
that relaxation vanishes in the Mean Field Approximation (MFA) and is slowed down as 
magnetization saturates when corrections to MFA are taken into account. 

For sake of simplicity we work with one of the most simple system which exhibits sym- 
metry breaking: the Ising model H = Jij^lSj. This model is related to the more 
familiar spin - boson model where spins interact with environment vibrations. If up and 
down local spin states are degenerated in energy, lattice degrees of freedom can be elim- 
inated by a canonical transformation. Ising model is then obtained as an approximation 
if lattice modes polarized in the spin quantization direction have smaller frequencies than 
perpendicular ones. Another physical realization is given by quantum dots interacting with 
the oscillation modes of the substrate. In this case the two states ( | t> ^"^^ I i> ) can 
be associated to the presence or absence of an extra electron. In Coulomb blockade condi- 
tion, i.e. when the two states are degenerated in energy, the coupling of quantum dots with 
phonons gives rise to an effective interaction between spins equivalent to the Ising model. 
Another possible physical realization of this model is a system of bosonic neutral atoms 
localized in the sites of a lattice interacting with repulsive interaction. Possible realization 
of such a system in an optical lattice has been discussed in 

We define the i-th qubit as the state \6i,(j)i >= cos(y)| ti> +6*"^' sin(y)| |j> associated 
to the corresponding site of a regular array ( i=0,..,N-l ). The local quantum coherence 
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operator is the projection operator on a given qubit state. 

7ri(^,0) = >< 0i,<^i| (1) 

The projection operator in terms of spin operators reads as 

TT,{9,(f)) = ^ + n,Si (2) 

where n is the unitary vector associated to the 9, (p direction and S is the vector spin 
operator. The statistical average of the quantum coherence operator is the qubit coherence 
probabihty (QCP) 

Pi is the probabihty of finding the spin at the i-th site, in the state 0j >, when the system 
is in equihbrium with an external bath. It is immediately seen from (|^) that as a consequence 
of the spin inversion symmetry of the Ising Hamiltonian H{{Si}) = H{{—Si}) we have pi = 
|. All directions 6i, (pi are then equally probable and coherence, as expected, is completely lost 
at equilibrium. Thus a "symmetric" thermodynamic state is obviusly useless as a reference 
state in a computation process. However, if a spontaneous symmetry breaking occurs below 
a critical temperature (Tc), we obtain a non trivial QCP. In our case spin inversion symmetry 
breakdown implies rrii =< >^ { angular brackets stand for thermodynamic average ). 
As a consequence the QCP becomes 

Pi = ^ + mi cos 9i (4) 

We see in that coherence is partially preserved at equilibrium even if the dependence 
on the phase (pi is lost. In the states space the direction 6*4 = is the most probable and 
the average of any even function of 9i does not vanish [^]. Hence, the thermodynamic 
state associated to the ordered spin phase is suitable as the reference state of a computing 
process. It is convenient to consider a reduced QCP as the difference between the ordered 
and disordered phase QCPs pi =< tTj >=< tTj — | >. We find convenient to identify the 
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0-logical state ( 1-logical state ) of the i-th site with the dependence of pi on cos 6*4 (sin^j). 
In other words different logical states are associated with different properties of a quantity 
related to the statistical average of the quantum coherence instead of different quantum 
states. 

Universal quantum computer needs 1-qubit and 2-qubits gates The first is a rotation 
of the 0-logic state in a superposition of and 1 logical states. For this reason we need a 
perturbation that modifies the local QCP introducing a term with depends on sin^j. The 
second modifies the two sites state introducing an entanglement between the two logical 
states. In the corresponding 2-qubits QCP a sin 6i sin 6j will appear as a result of the 
perturbation acting on a two sites logical initial state. 

For sake of simplicity we consider only impulsive perturbations. 

Ht = H + edit - Q+)Hi (5) 
After the action of a gate the QCP is given by 

Where 7fj(f) is the reduced quantum coherence operator in the Heisenberg representation. 

As 1-bit gate we consider a perturbation perpendicular to the z direction acting on a 
given site i=0. We choose to work with 

Hi = SS (7) 

If the Ising model is associated with localized particles ( quantum dots, optical lattices, etc. 
), the perturbation could be obtained by exchanging particles with a reservoir while in the 
spin case a magnetic field could be used. The new logical state is characterized as 

poii) = po cos e + Aoit) sin Oq sin e (8) 

where Ao^t) is given by 

Ao{t) = - sin ^0 < sin(2/iot + (/)o) > (9) 



and ho = J2j JojSj is the internal field. The strength of the perturbation determines the 
amplitude of the rotation from the 0-logical state (e = 0) and 1-logical state (e = |). Note 
that only the 1-logical state component of (|^) will evolve in time and then will be affected 
by decoherence { Ai{t) ^ ) while permanent modification are found in the 0-logical 
state amplitude. A simple oscillating behaviour is obtained in the MFA when internal field 
fluctuations are neglected, i.e. hi = Jiji^j and nii = ^ tanh It will be shown 

that a finite decoherence time arises as a consequence of non vanishing spin correlations 
Cij =< SiSj > ~ < Sj >< Sj > ( neglected in the MFA ). 

The 2-qubits gate can be realized by switching on a particles transfer between two sites. 
A simple result is obtained in the case of sites i,j belonging to two non-interacting adiacent 
layers of opposite magnetization. The Quantum coherence of the 2-qubits states associated 
with sites i and j is conveniently characterized by the statistical average of the product of 
reduced two sites quantum coherence operators pij =< (tTj — ^)(vrj — |) >. At equilibrium 
the two sites 0-logical states is characterized as 

Pij =< S^S] > cos Oi cos Oj (10) 

Note that in the MFA, where spin correlation are neglected, i.e. < SiSj >=< Sj >< Sj > 
{i 7^ j) the two sites reduced QCP is equal to the product of one site reduced QCP. The 
impulsive perturbation hamiltonian associated with sites i,j is then 

Hi = StSj + Sp- (11) 

Where i and j belong to two different non-interacting layers of opposite magnetization. In 
the time evolution of the reduced two sites QCP, following the impulsive perturbation we 
distinguish a constant term and a rotated one associated only with | \i.,[j> and | ii,tj> 
states. 

Pij (t) = Pij (0) + Aij sin (2e) sin 9i sin 6j (12) 

where Aij is given by 
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= j< (st - S]) sin [{hi - hj) t 



> 



(13) 



The impulsive local interaction produces the 1-logical state on two sites which oscillate with 
the interaction strenght. The peculiar superposition of and 1-logical states is a consequence 
of the quantum entanglement induced by the interaction. In the MFA limit, the internal 
fields hi and hj become non fluctuating quantities to be determined self consistently. To 
evaluate decoherence corrections to MFA must be taken into account in the framework of 
a perturbation expansion in the interaction range [1^. If the interaction range is small 
with respect to the system size the relaxation toward equilibrium can be seen as a quantum 
diffusion process wich spreads the perturbation over the whole lattice. We use the technique 
ref. of to calculate an auxiliary free energy functional where a single site time dependent 
field appears 



G"(A,t) = ln<^rr 



exp 



-(3H + Y.{2iatg, + Xk)S'^ 
k ) . 



(14) 



Here cr is a dicotomic quantity (cr = ±) and g], = Jq^ in the one site ( i=0 ) gate case and 
dk = {Jok — Jik) in the two sites ( i=0, j = l ) gate. In terms of this dynamical functional 
we can express both the one (P) and the two sites (|13D QCPs variations due to one and two 
sites gates 



^ = -^ aea;p(G"^(t)-G-(0)) 



cr=± 



(15) 



{A,}=0 



iij{(f>o-4>i) 



(dG" 


dG'^X 







{Aj}=0 



(16) 



As a simple generalization of procedure of ref. , we introduce a Legendre transform at 
fixed time and temperature 



(17) 



Where = The dynamical functional of ( P^ ) is obtained adding an imaginary time 
dependent term to the static functional. 
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F'^im, /3, t) = F{m, ^) 2zatg„ (18) 



k 



The time independent part of the (|T8D F{m,f]) is given in ref. |T0[- It is conveniently 
sphtted in two terms F{m,(3) = Fo(m) + Fi(m, /9) . The first term is associated to the 
non interacting system with fixed magnetization and the latter takes into account sites 
interaction. 



k 



{f^k + ^) ln(mfc + ^) + - ruk) ln(^ - m^) 



(19) 



Fi(m, (3) = +(3Y.Jik^i^k + E Jfk - ^fc') - "^'') + ■ ■ ■ (20) 

Ilk Ik 

The main point is to study the extremum condition |^ = in the complex plane. 

dF{m,(3) 



-2iatgk (21) 



For short time ( Jokt <^ 1 ) we can express the solution of the extremum problem as 
mlit) = nik + 5ml. where is the solution at t = 0. Expanding to the first order equation 
(pT|) we obtain 

^ omkOmi 



From ( ^21) we see that = itadnik where 5mfc oc 4 and, then, vanishes in the MFA limit 



The same expansion gives for the dynamical function of (|17| ) 

F"(m, (3, t) ^ F{m, f3) + J2 [2^^%"^fc - t'^Srrikgk] (23) 

k 

From (p2|) we can derive the time dependence of F{m) 

F^im, f3, t) ^ F{m, (3) + 2iatgkmk - Ft^ (24) 

k 

The quantity 2itaJ2k 9k''^k is related to the oscillating behaviour of the MFA while F = 
J2k Qk^T^k is real quantity related to the decoherence time. Because of the stability condition 
on the extremum of F{m), F > 0. If the number of next neighbors of site ( for one bit 



gate ), or sites and 1 ( for two bits gate ), which are next neighbors each other do not 
scale with the dimensionahty, then Srrik to first order in is equal to Srrik = fi'fc (| — 
and r becomes 

^-i:9l{\-ml) (25) 

As expected T is proportional to and vanishes in the limit of infinite dimensionality ( MFA 
limit ). It is important to observe that decoherence decreases as magnetization approaches 
saturation (m^ |). 

In conclusion we studied quantum computation in a spin system at thermal equilibrium 
in the presence of a spin symmetry spontaneous breaking. With respect to the conventional 
approach based on out of equilibrium systems two main results seem to be of interest. 

a) It is possible to define logical states as suitable thermal averages instead of quantum 
states. 

b) Decoherence, due to the conservative dynamics associated to energy operator wich 
defines statistical equilibrium, can be treated as a perturbation effect with respect to MFA 
and it is even possible to reduce it at will by increasing the order parameter. 

We think that the main results of our work should be valid for other systems which 
exhibit spontaneous symmetry breaking. 
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